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We study,numerically, the dynamical behavior of an electron in a two site nonlinear system driven 
by dc and ac electric field separately. We also study, numerically, the effect of electric field on single 
static impurity and antidimeric dynamical impurity in an infinite ID chain to find the strength of 
the impurities. Analytical arguments for this system have also been given. 
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I. INTRODUCTION 

The dynamics of an electron in presence of electric 
field, in a two level system as well as in a linear chain has 
been a subject of recent interest. (l]-|l| It is well known 
that for a double-well heterostructure in the absence of 
driving forces the electron can visit either well due to 
quantum tunneling. However, if the double well struc- 
ture is biased, i.e. driven by a dc field the electron in an 
eigenstate will be naturally localized in one of the wells 
because of the superposition of the coherent tunneling 
between two wells. || When the double well system is 
exposed to a pure ac field the electron gets frozen in the 
initially populated well, provided the ratio of the field 
strength to field frequency becomes a root of Bessel's 
function of order zero. || In presence of both ac and dc 
field, |6) the electron gets trapped in the initially popu- 
lated site when the ratio of the ac field strength to the 
frequency is such that the n th order Bessel function be- 
comes zero and the dc field strength becomes the n th 
multiple of the frequency of the ac field, provided the 
energy difference between the levels is much lesser than 
the frequency of the ac field i.e., J n (x n ) — 0,x n = E/uj, 

Eq = nco n = 0, 1, 2 and A/uj << 1, where A is the 

energy difference between the two levels. Here we con- 
sider the dynamics of an electron in a two site nonlinear 
system in presence of ac and dc electric fields separately. 

On the other hand a lot of studies have been done on 
the dynamics of a charged particle subjected to a periodic 
potential in presence of dc, ac as well as in the presence 
of both dc and ac electric fields. It is well known that, 
in presence of dc field, a charged particle executes oscil- 
latory motion in the reciprocal lattice, [jfj This, in turn, 
confines the particle spatially. || 

In presence of ac field, the particle is generally delo- 
calized except for the cases when the ratio of the field 
magnitude to the field frequency is a root of the ordinary 
Bessel function of order zero. |J In presence of both ac 
and dc field, the localization induced by the dc field can 
be suppressed by the ac field when the stark frequency 
becomes equal to an integral multiple of the frequency of 
the ac field. El 



Recently, Nazareno et. al [[TlJ had studied the dynam- 
ics of an electron in a ID chain in the presence of a single 
static impurity subjected to a dc electric field. They had 
shown that when the impurity potential ( eo ) is such 
that it coincides with the on-site energy due the field ( 
E ) on a particular site n, i.e., eEan — eo, where e is 
the charge of the particle and a is the lattice parame- 
ter, the packet oscillates resonantly between the impu- 
rity site and the particular site n. Therefore, when the 
electric field E — eo/ea, the charged particle oscillates 
resonantly between the impurity site and it's neighbor- 
ing site. However, we find that this condition is valid 
only if the particle is strictly confined to the impurity 
site and its neighboring site. On the other hand, in a 
ID chain there is a probability for the particle to tunnel 
to other sites. Because of this tunneling probability the 
resonance condition gets modified. 

We also study the dynamics of an electron in a ID per- 



fect chain in presence of two nonlinear impurities [|12 13 
which are next to each other, subjected to dc electric 
field. We suggest a possible means of measuring the 
strength of the nonlinear impurity. 

The organization of the paper is as follows. In Sec. II, 
we deal with a two site nonlinear system, in presence of 
electric fields. In Sec. Ill, we consider a one dimensional 
infinite chain with a static impurity and the dynamics of 
an electron in presence of a dc electric field. In Sec. IV, 
we consider a ID chain with two consecutive nonlinear 
impurity sites and the dynamics of an electron in presence 
of a dc electric field. Finally, in Sec. V we give a summary 
of our investigations. 



II. TWO SITE SYSTEM 

The time evolution of the electron in a two site nonlin- 
ear system in presence of both ac and dc fields is governed 
by the equations given by, 
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Here C\ and C 2 are the probability amplitudes of the 
particle to stay at site 1 and 2 respectively. V is the 
hopping amplitude of the particle between the two sites, 
Eq is the amplitude of the dc field, E is the amplitude 
of the ac field and u> is the frequency of the ac field. 
The kind of nonlinearity we consider here arises from the 
interaction of the electron with the vibration of the local 
oscillators in the system p2| . xi an d Xi are the nonlinear 
parameters determining the strength of the interaction 
of the particle with the local oscillators at site 1 and 2 
respectively. 

In the absence of both dc and ac electric fields (i.e. 
Eq = E = 0), the two site system is completely iso- 
lated from any external perturbations and the dynamics 
of the particle is governed completely by the strength of 
the nonlinearity at the two sites and the hopping matrix 
element between the two sites. 

In the special case when the strengths of nonlinearity 
at the two sites are equal (i.e. \i — X2 — x)> there ex- 
ists a critical value of x ( Xcr — 4V"), below which the 
particle remains delocalized between the two sites and 
above which the particle gets trapped at the initially pop- 
ulated site. [l2] | Thus, there is a abrupt transition from 
a completely delocalized state to a completely localized 
one across the critical value of the nonlinear parameter 
defined by Xcr- 

The other interesting situation arises when Eo = E = 
and xi = — X2 = X- This kind of nonlinearity appears 
when both the local oscillators oscillate in opposite phase 
and interact with the particle with similar strength. In 
this case, the time evolution of the probability of the 
particle at the initially populated site is given by, 
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In this limit, the transition is a continuous one, where, as 
the nonlinearity increases, the particle continuously gets 
trapped at the initially populated site. This is obvious 
from Eq. (J2J) where Pi(t) — > 1 as x - * 00 ■ 

We can, on the other hand, put xi = Xi — an d 
E = 0. This is the case when nonlinearity is absent 
and a pure d.c. field is applied to the two site problem. 
In this case, the two site problem effectively reduces to a 
two level problem, where the degeneracy of the two levels 
( due to equal site energies) is lifted by the externally 
applied dc electric field. It has been shown that, in this 
situation, the particle remains localized at the initially 
populated site provided the hopping matrix element is 
small compared to the strength of the electric field (i.e. 
V/E a «l). 

In the absence of nonlinearity and dc electric field 
(Xi = X2 = 0, Eq = 0) and in the presence of an ac 
electric field ( given by Ecos(uii)), the particle remains 
delocalized between the two sites except when the ratio 
of the amplitude of the a.c. filed (E) to its frequency (ui) 



is such that the zeroth order Bessel function Jq(E/u)) 
is identically equal to zero. Here too, the condition for 
localization is satisfied only when V/E << 1. 

The question naturally arises as to what happens to the 
dynamics of the electron in the nonlinear two site system 
in presence of field. Here, we consider the case where 
the nonlinear parameters are of equal strength but are 
opposite in sign. As mentioned earlier, this kind of situ- 
ation arises when the local oscillators oscillate in opposite 
phase but interact with the particle with same strength. 
Therefore, we consider xi = X an d X.2 — ~X- ln this sit- 
uation, the site energies of both the sites vary with time. 
Now, we consider the effect of a dc and ac electric field 
on the dynamics of an electron in the nonlinear two site 
system, separately. 

In presence of a pure dc electric field of strength Eq, 
the time evolution of the probability of the electron in 
the system is governed by the equations, 
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where P = pu- p 2 2, Q = P12 + P21, R = i{pi2~ P21) a nd 
Pij = CiC*, i,j=l,2. On solving Eqs. (||), the probability 
difference of the particle between the two sites is given 
by, 

PU P22 = (E ° n2 X)2 + ^COs(m) (4) 



where ft = y '4V 2 + (Eq — x) 2 - For the sake of discus- 
sion, we take \ an d Eo to be positive. Thus Eq. (Q) ex- 
plicitly shows that the particle becomes fully delocalized 
only when % = Eq. The field direction can be reversed 
by replacing Eq by —Eq. In this case however, there is 
no derealization for any field strength, since x IS still 
positive. 

We now consider the case of a pure ac electric field 
applied to the same nonlinear two site system. The dy- 
namics of the particle in this system is governed by the 
following equations. 

d -f = 2VR 
dt 



^ = -(Ecos(cjt)-x)R 

^ = -2VP+(Ecos(wt)- X )Q. 
dt 



(5) 



where P, Q and R are defined earlier. We cannot decouple 
these equations and hence solve them numerically, using 
fourth order Runge-Kutta method. We have checked the 
probability conservation of the particle at each time step. 
We find that, as long as V/E << 1, the particle becomes 
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fully delocalized only when x = nuj where n is a posi- 
tive integer and E/ui is such that the n th order Bessel 
function becomes nonzero (i.e. J n (E/u>) 7^ 0). On the 
other hand, the particle gets fully trapped at the initially 
populated site only when \ — nuJ an d E/u> is such that 
the n th order Bessel function is identically equal to zero 
(i.e. J n (E/cj) = 0). 




FIG. 1. The probability of the electron at site 1 (Pi) as a 
function of time, where, V = 0.8, ui = 2,E/uj = 1.9155 such 
that J\(E/oj) 7^ 0. The dotted curve, solid curve and the 
dashed curve are for ^=1.8, 2.0 and 2.1 respectively. 
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FIG. 2. The probability of the electron at site 1 (Pi) as a 
function of time, where, V = 0.8, a; = \ — 1,E/uj = 3.831, 
such that J\{E/uj) — 0. 



To show this, we have plotted P\ = -M^ ( probability of 
the particle to stay at the site 1) as a function of time in 
Figs. (0) and (§). In Fig. (§), we have taken V = 0.8, u = 
2 and E/u = 1.9155 such that Ji{E/uo) ^ 0. The dotted 



curve corresponds to % = 1.8 < w, the solid curve is for 
X = 2.0 = oj and the dashed curve for % = 2.1 > uu. 
The dotted and dashed curves show that, as long as the 
nonlinearity strength (x) is not an integral multiple of 
the frequency (u>) of the field, the electron remains mostly 
localized at the initially populated site. As soon as the 
nonlinearity parameter matches with the n th multiple of 
the frequency of the ac field (solid curve), the electron 
gets fully delocalized. On the other hand, in Fig. (Q), we 
have taken V = 0.8, ui = x = 1 an d E — 3.831 such 
that Ji(E/uj) = 0. The figure clearly shows that the 
electron gets completely frozen at the initially populated 
site. Thus, a two site system in presence of this kind of 
nonlinearity and ac field gives the same kind of dynamical 
behavior of the electron as is the case of a two level system 
subjected to dc and ac fields simultaneously. Here, the 
nonlinearity plays the role of a dc field in a two level 
system. 

The complete derealization of the electron in this sys- 
tem can be understood by the following physical rea- 
soning. Because of the presence of nonlinearity, the site 
energies of the two sites become different. Thus, the non- 
linearity maps the two site system to a two level system. 
Since site energy at both the sites depends on the prob- 
ability of the particle to stay at the respective sites, they 
fluctuate with time. However, the energy difference be- 
tween the sites, A = x(pn + P22) — X which is constant. 
This is because the total probability of the particle in 
the system is conserved. Therefore, when the particle 
is in one of the sites it can absorb an amount of energy 
X = nuj, where n is an integer, from the ac field and jump 
to the other site and thus, full derealization is obtained 
for this condition. 



III. ID CHAIN WITH SINGLE STATIC 
IMPURITY 

In this section, we consider a one dimensional infinite 
perfect chain with a single static impurity at the middle 
(0 th ) site. Our interest here is to find the dynamics of 
an electron initially populated at the impurity site in 
presence of a dc electric field. The governing equation 
for the electron dynamics in presence of the dc electric 
field is given by a set of differential equations, as, 



.dC n 
1 — — 
dt 



e n S n ,oC n + eaE nC n + V(C n+ i + C„_i), 

n = -N,-N + l, ,N-1,N. (6) 



C n is the probability amplitude of the particle to stay 
at the n th site, e„ is the site energy at the n th site, e is 
the electronic charge, a is the lattice spacing, E is the 
strength of the dc electric field and V is the hopping in- 
tegral between neighboring sites. N is chosen such that 
the particle does not feel the boundaries, which in turn 
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ensures the infiniteness of the chain. Without any loss 
of generality, we choose e = l,o = 1 and normalize all 
energies with respect to the hopping matrix element V. 
The delta function implies that the site energy of all sites 
except the zeroth site have been taken to be zero, in ab- 
sence of the dc electric field. To know the time evolution 
of the electron probability at site n we need to solve Eqs. 
(0). Since, it is not possible to solve them analytically, 
we use fourth order Runge Kutta method to solve them 
numerically. Here also, conservation of the probability of 
the particle has been checked at each time step. 
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FIG. 3. Time evolution of the probability of the electron in 
a ID chain with a static impurity at the th site, n represents 
the site index. The z-axis represents |C„| 2 . Here eo (strength 
of the impurity)= Eo (strength of the electric field) = 2. 



condition is larger than eo. We also find that, for smaller 
value of eo, the departure of this electric field strength 
(Eq ss ) from eo is larger and it decreases as eo increases. 
Thus, for very large values of eo, the resonance condition 
reduces to Eq cs — eo. We plot the time evolution of the 
probability of the electron in the system for Eq = eo = 2 
in Fig. (|^) and the same for Eq — 2.67, eo = 2 in Fig. 
(ky). It is obvious from Fig. (§) that for E = e = 2, the 
electron propagates through the system and gets con- 
fined among many sites. On the other hand, Fig. (Q) 
shows that the particle oscillates mainly between the im- 
purity site and one of its nearest neighbors. Therefore, 
for eo = 2, Eq — 2.67 is the electric field required for 
resonance between the impurity site and its neighbor. 
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FIG. 4. Same as Fig. (|) except E = 2.67. 

For a fixed value of e there is a particular value of the 
electric field strength {Eq cs ) at which the electron oscil- 
lates resonantly between the impurity site and one of its 
nearest neighbor sites. We observe that, when eo is small, 
then the electric field strength required for the resonance 



FIG. 5. Probability of the electron at site (Po) and 
at site 1 (Pi) as function of time. The solid curves are for 
eo = 4.0, Eo = 4.37 and the dotted curves are for eo = 4.0, 
E = 4.0. 

In Fig. (||) we have also plotted the probability of the 
electron at the impurity site and its nearest neighbor site 
as a function of time for eo = Eq = 4.0 (dashed curves) 
and e = 4.0, E = 4.37 (solid curves). Here too, it is 
clear, that for eo = 4.0 and Eq = 4.37 the particle os- 
cillates more strongly between the impurity site and its 
nearest neighbor than for eo = Eq = 4.0. Therefore, for 
eo = 4.0 one needs the field strength, Eq = 4.37, to get 
the resonant motion of the particle between the impurity 
site and its neighbor. In order to understand this de- 
viation from the normally expected resonance condition 
(i.e. eo = Eq) we find the resonance condition for general 
impurity strength as a function of electric field strength 
analytically. For this, we need to confine the particle be- 
tween the impurity site (site 0) and its nearest neighbor 
site (site 1). To do this, we renormalize the hopping ele- 
ments between the sites and -1, and between sites 1 and 
2 respectively. This in turn renormalizes the site energies 
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at sites and 1. The effective hopping between sites 
and -1 and that between sites 1 and 2 thereby decreases. 
The effective energies at site and 1 are given by 



4" = + V 2 /E Q 



E - V 2 /(2E ) 



(7) 



By putting e 
given by 



eff _ eff 



e-/ we get the condition for resonance 
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eo + y/4 + 6 



(8) 



This condition is in contrast to the condition for res- 
onance given by Nazareno et.al. |Tl|] Eq. (Q) explicitly 
shows that the resonance condition reduces to eo = Eq 
only when eo >> \/6. Thus, to get a resonant motion 
of the particle between the sites and 1, a higher value 
of electric field is required as compared to the impurity 
strength. We have plotted Eq £S — eo as a function of eo 
obtained from Eq. (||) (solid curve) along with the corre- 
sponding numerical plot ( dotted line) in Fig. (|J). 



FIG. 6. Eo es — eo as a function of eo- The dotted curve 
and the solid curve represent the numerical and analytical 
estimates, respectively. 



We see a good agreement between the numerical and an- 
alytical estimates. One can estimate the higher order 
corrections to Eq. (^) by further renormalizing the hop- 
ping matrix elements. However, the contribution from 
the higher order corrections is negligibly small. We there- 
fore, restrict ourselves only upto the first order correc- 
tion. Upto the first order, the strength of the static im- 
purity can be evaluated using Eq. (^]) , once there is reso- 
nant oscillation between the impurity site and its nearest 
neighbor, which is characterized by the emission of a pho- 
ton for the particular value of the dc electric field (E^ es ) : 
as remarked by Nazareno et.al. [[yj 



IV. NONLINEAR IMPURITIES IN A CHAIN 

In this section we consider a system of 1 dimensional 
chain with two nonlinear impurities at two consecutive 
sites. The nonlinearity arises due to the interaction of 
the particle with lattice vibration. We consider a special 
case where the particle interacts with the local oscillator 
only when it is at site or site 1. We denote the middle 
site as the th site. The interactions at both the sites ( 
site and 1) are equal in strength but opposite in nature. 
This may arise only when the oscillators at sites and 1 
oscillate in opposite phase. Assuming the presence of this 
kind of nonlinearity, we ask whether it is possible to find 
the strength of the interaction, by applying an external 
electric field, as was done for the single static impurity 
case in the previous section. Here, we give an answer to 
this question. 

We start by applying a pure dc electric field to the sys- 
tem and see the time evolution of the electron populated 
initially at the th site. The dynamics of the electron is 
governed by the equations given by, 

i^ = X \Co\ 2 C + V(C 1 + C- 1 ) 
at 

- (-x|Ci| 2 + So)d + V(C + C 2 ) 

at 



dC n 
i—— 
dt 



nE C n + V(C n+1 +C n ^), n^0,l. (9) 



We solve these equations numerically by using fourth or- 
der Runge Kutta method. We check the conservation of 
probability here too. We observe that, as E increases, 
the particle becomes more and more confined between 
sites and 1. The amplitude of oscillation between sites 
and 1 increases and then at some Eo ( say, Eq) it be- 
comes maximum. If Eq is increased further, the parti- 
cle gets trapped at the initially populated site. So, at 
Eq = Eq the particle oscillates between sites and 1 
resonantly. Here also, as is the case with a single static 
impurity in an infinite chain, for any strength of nonlin- 
earity, x the electric field needed to get the resonance 
condition is higher than \- This situation can also be 
detected by the radiation emitted by the electron while 
oscillating resonantly between sites and 1. 

We now find the resonance condition analytically. For 
this we have to reduce the system into an effectively two 
site system. In other words we need to confine the motion 
of the particle between sites and 1 . Since, at resonance 
condition, the time average probability of the particle at 
sites and 1 are approximately 1/2, we therefore replace 
the site energies arising from nonlinearity at site and 1 
by 7f and ^ respectively. We then renormalize the hop- 
ping matrix element between sites and its left neighbor 
and that between site 1 and its right neighbor ( the sites 
are numbered in increasing order from left to right, with 
the middle site being the zeroth site). As a consequence 
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of this renormalization, the site energies at sites and 1 
get renormalized and are given by, 




Now equating these effective site energies we get the the 
condition for the particle to oscillate resonantly between 
site and 1. This condition is given by 
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We plot Eq — x as a function of x obtained analytically 
in Fig. (Q) . Now we take a point from the curve for reso- 
nance condition (Fig. ^) which corresponds to x = 4.0 
and Eo = 4.34. For these values of x an( i Eq we see the 
the time evolution of the electron probability in the sys- 
tem. This is shown in Fig. (||). We notice that the elec- 
tron mostly oscillates between sites and 1 resonantly. 
This also confirms our analytical estimate for resonance 
condition. Thus tuning the electric field one can obtain 
the resonance condition and hence the strength of the 
nonlinear impurity can be estimated using Eq. ( [TT| ) pro- 
vided that special kind of nonlinearity is present in the 
system. 

V. CONCLUSION 

We have studied the dynamics of an electron in a two 
site nonlinear system in presence of externally applied 
dc and ac electric fields, separately. We have shown the 
localization-delocalization conditions in both the cases. 
We have also studied the dynamics of an electron in a in- 
finite ID chain having a single static impurity in presence 
of a dc electric field. We have given an analytic expression 
for the condition when the particle oscillates resonantly 
between the impurity site and its nearest neighbor. This, 
in turn, gives a possible method of measuring the impu- 
rity strength. A similar analysis has been performed for 
an antidimeric dynamical impurity present in an other- 
wise perfect chain. 



FIG. 7. Eq — x as a function of %■ 
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FIG. 8. Time evolution of the probability of the electron in 
a ID chain with nonlinear impurity \ at site and -\ at site 
1. n represents the site index. The z-axis represents |Cn| 2 . 
Here X = 4.0 and E = 4.34. 
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